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Abstract 

Polarimetry of the Cosmic Microwave Background (CMB) represents one of the possible di- 
agnostics aimed at testing large-scale magnetism at the epoch of the photon decoupling. The 
propagation of electromagnetic disturbances in a magnetized plasma leads naturally to a B-mode 
polarization whose angular power spectrum is hereby computed both analytically and numer- 
ically. Combined analyses of all the publicly available data on the B-mode polarization are 
presented, for the first time, in the light of the magnetized ACDM scenario. Novel constraints on 
pre-equality magnetism are also derived in view of the current and expected sensitivities to the 
B-mode polarization. 
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1 Formulation of the problem 



When (linearly) polarized electromagnetic radiation propagates in a diamagnetic material, the polar- 
ization plane of the wave changes by an angle A$(A, B\\) 



where V denotes the Verdet constant, B\i characterizes the magnetic field intensity along the direction 
of propagation of the electromagnetic radiation of wavelength A; L is the distance travelled by the 
wave within the medium. The Verdet constant V not only depends upon the wavelength but also on 
the specifics of the given material such as, for instance, the temperature and the charge concentration. 
In electromagnetic plasmas, the kinetic energy of the charge carriers dominates against the potential 
energy, and the Verdet constant can be computed in simple terms: 



where mg is the electron mass, fie is the charge concentration of the electrons and the integral allows 
for the spatial variation of the magnetic field intensity as well as of the charge concentration. Equation 
( 11.21) holds under two physical assumptions (see, e.g. [Ill2]): 

• the plasma is globally neutral, i.e. the charge concentrations of the electrons and of the ions are 
balanced, i.e. fie = fii = ng (where uq denotes generically the common value of he and hi); 

• the plasma is cold, i.e. the kinetic temperature of the electrons and of the ions is always much 
smaller than the corresponding masses. 

Since the plasma is cold (i.e. non-relativistic) and globally neutral the values of the plasma and Larmor 
frequencies for the electrons are always smaller than the corresponding quantities but computed for 
the ions: this is why, in Eq. (11.21) only the electron mass and the electron concentration appears. 
In the case of a tokamak Eq. (11. 2p is often used to determine the charge concentration he- Linearly 
polarized electromagnetic radiation can be sent through the tokamak and the A$(A^) can be directly 
measured, for instance, as a function of the frequency of the incident radiation. The latter experimental 
information is sufficient to determine, at least approximately, the electron concentration he since the 
geometry and intensity of the magnetic field are known from the design of the instrument. 

Given the profile of the electron concentration it is relevant, in diverse circumstances, to assess 
the magnetic field intensity. This is, in a nutshell, one of the objectives of cosmic polarimetry: the 
polarization properties of the cosmic microwave background radiation (CMB in what follows) are 
exploited as a diagnostic for the existence of large-scale magnetic fields. The polarization observable 
of the radiation field are customarily parametrized in terms of the Stokes parameters [3] whose evolution 
can be written in a schematic form which is closely analog to the set of equations often employed to 
describe the propagation of polarized radiation through stellar atmospheres [H E] : 









dr 




(1.4) 



(1.5) 



where 
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• e' is the differential optical depth of the system (mainly due, in the CMB case, to photon-electron 
scattering) and K-p is the Faraday rotation rate, i.e. /tp = 3 (n ■ -B)/(16 vr^ e); 

• Si and Sq are the source terms for the intensity and for the polarization; 

• within the present notations, the (rotationally invariant) / and V Stokes parameters describe, 
respectively, the intensity of the radiation field and the circular polarization; the Q and U Stokes 
parameters are not rotationally invariant and are both sensitive to the linear polarization. 

Furthermore, according to the standard conventions, we also have: 

1 fU 



$ = -arctan(^-j, P = ^Q^ + U^. (1.6) 

Equations f ll.3p - fll.5p are the flat-space analog of the evolution equations studied in the present paper. 
In Eqs. (11.31) . (11.41) and (11.51) the convective derivative with respect to r is actually a total derivative 
with respect to the cosmic time coordinate which receives contribution also from the fluctuations of 
the geometry. 

The large-scale fluctuations of the geometry induce computable inhomogeneities on the Stokes pa- 
rameters and they will be denoted, respectively, a^ Ai(n, r), AQ(n, r) and Au(n, r). The experimental 
observations on the properties of the CMB polarization are expressed not in terms of Aq (n,r) and 
Au (n,r) but rather in terms of the corresponding E-mode and B-mode polarizations. The rationale 
for this common practice stems directly from the observation that the brightness perturbations of the 
Stokes parameters sensitive to linear polarization (i.e. AQ(ri,r) and Au(5t., t)) are not invariant under 
rotations on a plane orthogonal n. Consequently, their orthogonal combinations 

A±(n,r) = AQ(n,r)±2Au(n,r), (1.7) 

transform, respectively, as fluctuations of spin- weight ±2 [7j (see also [8j). Owing to this observation, 
A±{n,T) can be expanded in terms of spin-±2 spherical harmonics ±2Yem{n), i.e. 

A±(n, t) = Y^ a±2,im ±2Yim{n). (1.8) 

em 

The E- and B- modes are, up to a sign, the real and the imaginary parts of a±2/m, i-e. 

In real space (as opposed to Fourier space), the fluctuations constructed from and a}J^ have the 
property of being invariant under rotations on a plane orthogonal to n. They can therefore be expanded 
in terms of (ordinary) spherical harmonics: 

AEin,T)=J2Ni'a2YUn), AB(n, r) = J] iV^^ aS r,„(n), (1.10) 

Im tm 



where = ^{i - 2)\/{i + 2)1. Under parity inversion aS ^ (-l)^aS while af^ ^ (-l)^^^aS 
The EE and BB angular power spectra are then deflned as 



^In the spherical basis h = (i?, ip) and dh = sind M dip = —dcos.ddip. The notation /i = cosd wiU be adopted. 
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where (...) denotes the ensemble average. In the minimal version of the ACDM paradigm the adiabatic 
fluctuations of the scalar curvature lead to a polarization which is characterized exactly by the condition 
0^2, im = 0,-2, im, i-G. a)„ = 0. In more general situations also the cross-correlations of the various modes 
can be defined in full terms. Having introduced the appropriate rotationally invariant observables, the 
empirical observations can now be quickly summarized, in a buUeted form, aj^: 

• the TE and the EE angular power spectra have been measured so far by various experiments 
with different accuracies; 

• direct measurements of the BB angular power spectra led just upper limits on the B-mode 
polarization since the observations are consistent with the null result; 

• the polarization power spectra detected by different experiments are consistent. 



Experiments 


Years 


Data 


Pivot frequencies 


Upper limits on BB 


Dasi la m [lU 


2000-2003 


EE, TE, BB 


26-36 GHz 


<50 (/iK)2 


Boomerang [12] 


2003 


EE, EB, BB 


145 GHz 


e^f"^ <8.6(/iK)2 


Maxipol [ni [m 


2003 


EE, EB, BB 


140 GHz 


^f^^ < 112.3 (/iK)2 


Quad [151 [m [m 


2007 


EE,TE, TB 


100/ 150 GHz 




cbi [n] 


2002-2005 


EE, TE, BB 


26-36 GHz 


^f^^ < 3.76 (^K)2 


Capmap [19] 


2004-2005 


EE, BB 


35-46/84-100 GHz 


e^f"^ <4.8(/iK)2 


Wmap [201 EI] 


2001-2006 


EE, TE, BB 


23-94 GHz 


g?""^ < 0.15 (/iK)2 



Table 1: The main polarization experiments, their typical frequencies, their salient observations are 
illustrated in the light of the polarization properties of the CMB. In the last column Q^^^^ = £{£ + 
l)CrV(2vr). 



In Tab. [T]the salient informations on different experiments measuring the polarization properties of the 
CMB are swiftly collected. In Tab. [H the last two columns at the right include specific informations 
on the frequency channels characteristic of each experiment and of the upper limits on the B-mode 
autocorrelation. The frequency of the observational channel constitutes a relevant piece of information 
since, as already remarked, the Faraday rate depends upon the wavelength of the propagating radiation. 
The main references to each experiment have been collected in the first column of Tab. [TJ In the case 
of the WMAP experiment the most recent data release can be found in [211 1221 [231 1211 [25] . The first 
data release corresponds to [Sni [261 EZ] • For the specific aspects of the polarization measurements it 
is also interesting to bear in mind Ref. [2S] (see also 

The paradigm in vogue these days is the ACDM scenario where A stands for the dark energy 
component and CDM for the cold dark matter component. In the ACDM context the origin of the linear 
polarization of the CMB stems directly from the properties of last scattering in conjunction with the 
presence of a primordial spectrum of curvature perturbations which is said to be adiabatic. Adiabatic 
means, in the present context, that the curvature perturbations do not arise as inhomogeneities in 

^Following a common practice the polarization will be described in terms of the autocorrelations and cross-correlations 
of the E-modes and of the B-modes (i.e., in the jargon, the EE, BB and EB angular power spectra). Relevant pieces 
of information are also encoded in the cross-correlation between the temperature and the polarization, i.e. the TE and 
TB correlations. 
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the plasma sound speed but, on the contrary, they are related to genuine fluctuations of the energy 
density and of the spatial curvature. In the adiabatic case the location of the first acoustic peak of 
the TT angular power spectrum determines the location of the first anticorrelation peak of the TE 
spectrum. This is what is observed since, roughly, £tt — 220 and £te — 3£tt/4 where £tt denotes 
the location of the first acoustic peak and £te corresponds to the first anticorrelation peak in the TE 
spectrum. Finally the success of the ACDM scenario in pinning down the observed features in the TT 
and EE angular power spectra gives us also, indirectly, some precious informations on the dynamics 
of photon decoupling and on the possibility of late-time reionization. As indicated by Tab. [1] various 
experiments measured (and are now measuring) the EE and TE angular power spectra. On top of the 
5-yr WMAP results [211 [22l [231 1211 [25] recently there have been rather promising results coming from 
the Quad collaboration [16] (see also [151 [H]) claiming rather intriguing measurements of the EE and 
TE angular power spectra up to rather large harmonic^, up to £ ~ 1500 and even i ~ 2000. 

In the ACDM scenario with no tensors, the B-modes vanish except for a minute component coming 
from the gravitational lensing of the primary anisotropies. If the ACDM scenario is complemented by 
means of a tensor contribution, then the value of the BB autocorrelations will depend upon the ratio 
between the power spectrum of the tensor fluctuations and the power spectrum of the fluctuations of 
the scalar curvature. The current experimental data reported in Tab. [1] are also illustrated in Fig. [H 
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Figure 1: Measurements of the BB autocorrelations as they have been reported in the various experi- 
ments listed in Tab. [H 

While the possible presence of a B-mode polarization can be attributed to a stochastic background of 
relic gravitons, Eqs. (I1.3p - (I1.5I) suggest that it could be the signal of a pre-decoupling magnetic field. 

Equation (11.31) describes the evolution of the intensity of the radiation field and, in a magnetized 
medium, the source term Sj receives contribution from both the fluctuations of the geometry as well 
as from the fluctuations of the plasma which is, in turn, sensitive to the stochastic magnetic field itself. 
The remaining two equations determine the polarization observables and can be solved by iteration, 

^The given measured multipole £ is related to the angular separation. Large multipoles correspond to small angular 
separations, i.e., roughly d — n/i. 
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i.e. by positing that 

g = g + g«, t/ = F + f/«. (1.12) 

The first iteration determines Q which is a solution of Eqs. fll.3l) - fll.4l) with [7 = 0, i.e. 

^ + e'g = e'^Q, ^ + e'f/(i) = -2e'«:Fg. (1.13) 

The result of the second iteration produces U^^^ ^ and, ultimately, a B-mode polarization. In the 
present paper the magnetized B-mode will be computed both analytically and numerically. Further- 
more various bounds on the magnetic field intensity will be derived by comparing the Faraday-induced 
B-mode both with the upper limits on the BB angular power spectra (see Tab. [T]) as well as with the 
pubhshed data on the B-mode angular power spectrum (see e. g. Fig. [T]). 

The present investigation is based on a number of previous results. The large-scale magnetic fields 
are included both at the level of the evolution equations as well as at the level of the initial conditions. 
This program has been undertaken in [30] (see also [31] |32j) and subsequently developed in [331 El]- 
Semi-analytic results for the polarization observables which are relevant for the present investigation 
have been discussed in [33]. Recently the results of [30l[33l[34] have been complemented by a dedicated 
numerical approach [3S1 ESI EZ] which is based, indirectly, on COSMICS [5HI EH] and, more directly, 
on CMBFAST [IQIIIT]. 

In [27] a bound on pre-decoupling magnetic fields has been derived, by including, as suggested in 
[33], the large-scale magnetic fields in the heat transfer equations as well as in the initial conditions 
of the full Einstein-Boltzmann hierarchy. The obtained results showed that the bounds obtainable 
from the TT correlations were still more constraining than the direct Faraday bounds stemming from 
the upper limits on the B-mode autocorrelations. The assumptions of [2Z] will be relaxed and the 
finite thickness of the last scattering surface will be taken into account in the analytical and numerical 
estimate of the Faraday autocorrelation. While the sudden decoupling limit (i.e. infinitely thin last- 
scattering surface) is justified for large angular scales (i.e. i < 100) over small angular scales finite 
thickness effects are comparable with diffusive effects. The latter statements can be understood also 
in analogy with the calculations of the temperature autocorrelations. While the finite thickness effects 
can be ignored for the estimate of the Sachs- Wolfe plateau, the same cannot be said for the region of 
the acoustic peaks where both the finite thickness effects compete with diffusive damping to determine 
the well known shape of the temperature autocorrelations. Semi-analytic treatments of the finite 
thickness effects can be found in [53l [55l [5ll [561 EI]- Finally, finite thickness effects also determine 
the correlated distortions of the first three acoustic peaks in the temperature autocorrelations when 
large-scale magnetic fields are consistently included in the calculation [311 ES] • 

In the light of some of the recent (and forthcoming) experimental data, finite thickness effects play 
a relevant role in the estimate of the Faraday autocorrelations as well as in the estimate of the B-mode 
autocorrelations. Indeed (see Tab. [Tjand Fig. [T|) the upper limits on the B-mode autocorrelations stem 
from rather small angular scales. Over those scales the Faraday induced B-mode cannot be trusted in 
the sudden decoupling limit. 

While physically the small angular scales are the realm of diffusive effects, geometrically the small 
angular separations imply that the microwave sky degenerates into a plane. On a 2-sphere the angular 
eigenfunctions of the Laplacian are spherical harmonics. When the radius of the sphere is taken 
to infinity (i.e. the angular separation becomes small) the eigenfunctions of the Laplacian become 
plane waves. In the language of group representations this is one of the simplest examples of group 
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contraction (in the Inonii-Wigner sense) where the group SO (3) contracts to M(2), i.e. the semidirect 
product of the Euchdean group in two dimension and of the SO (2) group [58]- This observation can 
simphfy the evaluation of comphcated correlators and will lead, ultimately, to an analytical expression 
of the Faraday-induced autocorrelations at small angular scales. 

The plan of the present investigation is the following. In section 2 the essentials of the interplay 
between large-scale magnetism and CMB physics will be reviewed. Faraday autocorrelations will be 
computed in section 3. The B-mode autocorrelations will computed in section 4. In Section 5 the 
parameters of the magnetized ACDM model will be analyzed within a frequentist approach. The 
concluding remarks are collected in section 6. 



2 Pre-decoupling plasma and large-scale magnetic fields 

Heeding observations, the linear polarization of the CMB is a well established empirical notion. Con- 
versely the optical activity of the electron-baryon plasma is still under active scrutiny [32] • The 
program suggested in [32] has been to include consistently large-scale magnetic fields in all the steps 
leading to the angular power spectra of the CMB anisotropics. One of the hopes of the program will 
be to gauge accurately the effects of large-scale magnetic fields and to set limits on their putative 
existence. The essentials of the theoretical framework will now be schematically introduced. More 
detailed derivations can be found in [30l [331 1311 [351 EZJ- The present approach to pre-decoupling 
magnetism is rooted in three plausible assumptions: 

• [i] large-scale magnetic fields, if present, do not break the spatial isotropy and, therefore, cannot 
be uniform; 

• [ii] the curved-space effects as well as the relativistic fluctuations of the geometry are described 
within the lore provided by General Relativity; 

• [iii] the same physical laws describing the dynamics of electromagnetic fields in terrestrial plasmas 
are fully valid. 

The implications of each of the three aforementioned hypotheses will now be briefly introduced with 
particular attention to the tools which will be exploited in the subsequent part of the present analysis. 
The necessity of leaving unbroken spatial isotropy demands that the Fourier modes of the magnetic 
field must obey 

{B,{k,T)B,{p,T)) = — VB{m,{k)6^'\k + p), (2.1) 

where 

P^^ik) = (5,, - PB(fc) = ^b(|-)""'. (2.2) 

In Eq. (12.21) As is the spectral amplitude evaluated at the magnetic pivot scale /cl; is the magnetic 
spectral index. Within the conventions of Eqs. ( 12. ip and (12. 2p . As will have dimensions of an energy 
density and will therefore be measured in nC^. The scale- invariant spectrum of magnetic inhomo- 
geneities does correspond to ne = 1. In Eqs. (I2.1I) - (I2.2I) the correlation functions have been assigned 
within the same conventions employed to assign the curvature perturbations on comoving orthogonal 
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hypersurfaces, i.e. 

{n.ik)n.ip)) = j^Vnik)5^'\k + p), Vn = An[-j . (2.3) 

In Eq. ( 12. 3p A-ji denotes the amplitude of curvature perturbations at the pivot scale kp = 0.002 Mpc^^. 
Owing to the isotropy (and spatial flatness) of the background geometry, the line element can be 
written, in the conformally flat parametrization, as 

ds^ = a^{T)[dT^ -dx% (2.4) 

The scale factor will always be normalized in such a way that oq = 1. With the latter convention 
comoving and physical scales coincide at the present time. This convention will be enforced both in 
the analytic calculations as well as in the numerical applications. 

Using Eq. (12. 4p . the second hypothesis listed above (i.e. [ii]) implies that a(r) obeys the Friedmann- 
Lemaitre equations which can be written as 

-n' = 47iGa\p, + pt), (2.5) 
a' 

0, H=-, (2.6) 
a 

where pt and pt denote, respectively, the total pressure and the total energy density, i.e. 

Pt = Pe + Pi + Pu + + Pc + Pa- (2.7) 

The various components of Eq. fl2.7p refer, respectively, to the electrons, the ions, the neutrinos, the 
photons, the cold dark matter (CDM) species, and the dark energy. The relativistic fluctuations of the 
geometry will be described within the longitudinal gauge where the perturbed entries of the metric 
take the form: 

Ssgoo = 2a^4){x, r), d^gij = 2a^i){x, T)6ij. (2.8) 

As previously stressed (see, e.g. [301 [331 [36]) the longitudinal gauge is particularly suitable to carry on 
semi-analytic estimates. Conversely the synchronous coordinate system turns out to be rather handy 
in dedicated numerical calculations. The subtleties related to the use of different coordinate systems 
in the description of magnetized plasmas have been thoroughly explored in the recent past pT] (see 
also [36]). 

From the third (and last) assumption listed at the beginning of this section the evolution of the 
electromagnetic fields follows, schematically. Maxwell's equations 

V ■ E = 4:7re J dh[f+{x,v,T) - f_{x,v,T)], V ■ 5 = 0, (2.9) 
VxE + B' = 0, V X B- E' = 4:TTeJ d^vv[U{x,v,T)- f^{x,v,T)], (2.10) 

where f±{x, v, r) are the one-body distribution functions for electrons and ions written in terms of the 
comoving velocity v (and not upon the comoving three-momentum q) since, after neutrino decoupling, 
electrons and ions are non-relativistic, i.e. rrie ^ and mj ^ Tj where Tg and denote, respectively, 
the electron and ion temperatures which are, in the first approximation, of the order of the temperature 



2 SttG 2 



-a Pt, 



p[ + m{pt + pt) 
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of the photons. Consequently the comoving three-velocity and the comoving three-momentum are 
related as 

Ve,i = I '^"^ ^ (2.11) 

^jql^ + ml^a? "^eiO- 

where the subscripts refer, indifferently either to electrons or to ions. In Eq. (12.111) the relativistic 
fluctuations of the geometry have been neglected just to avoid complications. In explicit calculations, 
however, they must be taken into account. Indeed, the one-body distribution functions will obey the 
corresponding Vlasov-Landau equations, i.e. 

^ + ± + ^^B.e^"')^-^ + - ""^^^^^ = (2.12) 

where ip and are the two entries of the perturbed metric introduced in Eq. (12.81) and Ccou is the 
collision term which is dictated by electron-ion scattering as well as by the scattering with the photons 
whose concentration, around decoupling, is lO^'' times larger than the one of electrons and ions [37j . 
The electromagnetic fields appearing in Eqs. (I2.9I) - (I2.10I) are all comoving but, as it can be appreciated 
from Eq. (12.111) the presence of electrons and ions breaks explicitly the conformal invariance of the 
whole system. In particular this will reflect, as discussed in the past, in the nature of the tight-coupling 
approximation [3T] . 

Equations (12.91) . (I2.iup and (12.121) summarize, in a nutshell, some of the salient features of the 
description of electromagnetic plasmas in curved space-times and in the presence of the relativis- 
tic fluctuations of the geometry. At early times (i.e. well prior to matter-radiation equality) the 
Boltzmann-hierarchy must be treated by making explicit the tight Coulomb coupling as well as the 
tight Thomson coupling. Later on this approximation will break down. The general approximation 
scheme will now be outlined leaving aside details which can be however found elsewhere. 

Given the minuteness of the Debye scal^ in comparison with the Hubble radius the plasma ap- 
proximation works extremely well. A two fluid description can then be employed to describe the 
electron-ion system before decoupling. Taking the first moment of Eq. (12.121) we can get the evolution 
of the comoving velocity fields for electrons (i.e. Ve) and for ions (i.e. Vi): 

v,' + Hv, = '-[E + v,xB]-V(P + Cei + Ce^, (2.13) 



Vi' + Hvi = — j[E + ViX B]-V(f) + Cie + Ci^, (2.14) 



where = a^Ug and rii = a^hi. In Eq. (I2.13P Cd and Ce-y denote, respectively, the collision terms 
between electrons and ions and the collision terms between electrons and photons. With analog 
notation the collision terms of the ions are included in Eq. (12.141) . The electron-photon cross section 
dominates the mean free path of the photons because of the hierarchy between the electron and proton 
masses. 

For comoving scales much larger than the (comoving) Debye scale Ad and for frequencies much 
smaller than the plasma frequency the electromagnetic fields obey an effective one-fluid description 
which arises naturally by assuming that the Coulomb coupling between electrons and ions is strong 



^The Debye scale has, in the present context, exactly the same meaning and the same form it would have, for instance, 
in a terrestrial plasma, i.e. Ad — \/Tc/{8Tr e'^he). For typical lenght-scales L > Ad the collective effects dominate and 
the positive charge carriers are compensated by the negative charge carriers. In the opposite limit the propagation of 
electromagnetic disturbances cannot be neglected. 
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. In this limit the effective set of equations is given by the generahzation of magnetohydrodynamics 
(MHD) to the case when the geometry is curved and it fluctuates according to Eq. (12.81) . This set 
of equations has been written in [33] and will not be repeated here. The evolution equation of the 
electron- ion velocity can be directly obtained by carefully summing up Eqs. fl2.13p and (12.141) . After 
simple algebra, the result ij^: 

e'^ + = _ + ^Pie'ie, - e^), (2.15) 

where 6'^ = V ■ "iL is the three-divergence of the comoving three- velocity of the photons and wher^: 



e = rieXeCTTha, cTTh = ^vr( — ) , (2.16) 

Oh = : , 6'e = V-We, 6'i = V-Wi. (2.17) 

rric + rrii 

In Eq. (12.151) . J denotes the (solenoidal) Ohmic current. Notice that is effectively the centre-of mass 
velocity of the electron-ion system; 6^ in the limit of strong Coulomb coupling 6^ coincides with the 
bulk velocity of the plasma. Equation (I2.15P is coupled directly both to the evolution of the density 
contrasts as well as to the evolution of the photon velocity field, i.e. 

S'^ + p^-iij' = 0, (2.18) 

S'^ + Oh- 3^' = f "'^ . , (2.19) 



0; + ^V25^ + VV = e'(^^b-^7). (2.20) 



Since the electron-ion system is strongly coupled also to photons it is in principle possible to consider 
photons, baryons and electrons as a unique dynamical entity [33] : 



+ "l^".^ - 4,..^ (1 ^ J,,) V ■ X i.) - - (2.21) 

wher^ Rhir) 

R^(^r) = -P^ = 690.18 f-^^^V(^), (2.22) 

^ ' Ap^ V0.02273/ ^ ^ ' 

is the baryon-to-photon ratio. Equation (12.211) can be easily derived by combining Eqs. (12.151) and 
(I2.20p and by noticing that when the photon-electron-baryon coupling is tight 6^ — 0^ c:^ 9^^. By 
going from Eq. (12.120 to Eq. (12.220 the effective number of degrees of freedom has been progressively 
reduced. The effective description mentioned so far is very useful at early times and over large scales 
(i.e. small multipoles). It should not be applied, however, outside its domains of validity. 



^Typically toi ~ m-p where m-p is the proton mass. 

^The absence of the ions from the photon mean free path (and in the corresponding rate) is due to the fact that the 
ions contribute httle to the photon mean free path since their cross section is neghgible in comparison with the one of 
the electrons. 

^FoUowing the usual notations ilbo = Pw/Pcrit denotes the critical fraction of baryons at the present time. For a 
generic component x we will have that Q^q = Pxo/ pmt- 
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The breakdown of the effective one-fluid description is related to the possible propagation of elec- 
tromagnetic disturbances. Indeed, in MHD the Ohmic current is solenoidal (i.e. AttJ = V x 5) and 
the electric fields vanish in the baryon rest frame because of the large value of the conductivity [30j. 
The latter statement depends upon the frequency of the electromagnetic radiation. It is actually well 
known that high-frequency waves can propagate inside (magnetized) globally neutral plasmas [HE]- 

The baryon-photon system ceases to be a unique physical entity around the decoupling time. As 
time goes by the velocity fields of the photons and of the baryons must be separately followed by 
numerical means. The electron-ion system does not break down after a particular time but above 
a particular frequency of the spectrum of plasma excitations. If the plasma is not magnetized the 
only typical scale of the problem is set by the electron and ion concentrations which determine the 
corresponding (comoving) plasma frequencies, i.e. 



^pe = \ = t^pefl, ^pc = \ , (2.23) 

V rripa \ 



AirrioXee^ _ 

pi — \ = ^pi'^) ^pi \ 

\ mpa \ 



«.a;.e^^ (2.24) 



where it has been used that no = a^Ue = a^Ui. Taking into account that the comoving concentration 
of charge carriers is 10^" times smaller than the one of photons 

"» = ^.".o. -)« = 6.219 X W-^i^) {^)\ (2.25) 

we can say that the effective bulk description provided by Eq. (12.151) and its descendants break down 
for physical frequencies larger than [3] 



u;pe(z) = 0.0285 V^f-^^^^^y^Yi±l') MHz, Jz) = Up,{z) . (2.26) 

' ^ VO.02273/ VlOOOy ' V "^p 

As already stressed in [37] to treat the dispersive propagation of electromagnetic waves in a plasma 
demands necessarily that electrons are treated separately from the ions. The second point to be 
emphasized is that the inhomogeneity scale of the magnetic fields is typically much larger than the 
Larmor radius. This hierarchy of scales justifies the use of the guiding centre approximation [U [2] in 
the derivation of the relevant dispersion relations. The guiding centre approximation is well justified 
since, in the problem at hand, the scale of variation of the magnetic field is (always) much larger than 
the gyration radius for ions and electrons. 

The spectrum of plasma excitations may be in principle very rich and it includes a load of different 
phenomena which have both practical and conceptual implications |42]. The dispersive propagation 
associated with the ordinary and extraordinary wave is, however, irrelevant in the case of the CMB 
|37j . In similar terms, other dispersive effects are also absent. The rationale for these occurrences 
stems from the hierarchy between the plasma (and Larmor) frequencies of electrons and ions (see e.g. 
Eq. fl2.26p ) and the frequency of the maximum of the CMB which is given by 

^^max = ^ T^o = 222.617( ^^^°^ ) GHz, w = 3.920, (2.27) 

where x = k/T^o and Xmax corresponds to the maximum of the spectral energy density of the CMB, 
i.e. ^ 

o (x) = l^o o^L^, o = 2.47 X 10"^ /iq I (2.28) 



10 



By taking the derivative of Q^{x) with respect to x, the maximum of Q^{x) is determined by the 
condition [4 + e^{x — 4)] = which imphes that Xmax = 3.920. The angular frequency related to i^max 
is given by cUmax = 27rz/max- By comparing cUmax with the other typical frequencies of the problem we 
have that 

Wmax > tUpo > t^Bc > ^^pi > t^Bi, (2.29) 

possible resonances in the dispersion relations are then avoided, but, still the polarization plane of the 
CMB can be Faraday rotated by an angle 

To) = A{V) r /C(r) \h ■ B{x, t)] dr, /C(r) = e'(r)e~^("'"«). (2.30) 
Jo 

In Eq. (12.301) the squared brackets have been used to remind of the vector structure of the integrand; 
A{T') = 3l7~^/(167r^e) depends upon the comoving frequency V; finally in Eq. (I2.30p /C(t) is the 
visibility function expressed in terms of the optical depth e(r, tq) 

/TO 
Xc{r') CTTh nc(r')a(r') dr'. (2.31) 

The visibility function vanishes for r ^ Tj-cc and has a maximum around recombination, i.e. when 

e" + e'2 = 0, ^ = -a;e(r)ne(r)aTha(r) = -t' . (2.32) 
dr 

The second expression in Eq. (12.321) clarifies that a minus sign appears in the time derivative of e(r, tq) 
since r appears in the lower limit of integration. The magnetic fields not only enter the evolution 
of the charged species and of the photons (which are tightly coupled to electrons and baryons before 
decoupling) but also the (perturbed) Einstein equations. The Hamiltonian and momentum constraints, 
stemming from the (00) and (Oi) components of the perturbed Einstein equations are: 

V'V' - + ^') = A'nGa^^pt + <5sPb], 4pB(r,f) = /'^f (2.33) 



(Pt + Pt)Ot + ^ ' 



(2.34) 



47ra4 

In Eqs. (I2.33P and fl2.34p d^Pt and 6't, denote, respectively, the total density fluctuation of the fluid 
mixture and the divergence of the total velocity field (i.e. 9^ = div\) whose expressions, in terms of 
the four components of the plasma, i.e. i/, 7, c (CDM) and b (baryons), is 

(^sPt = X! '^sPa, 4Pt = SsPa., {Pt + Pt)Ot = XI (Pa + Pa)6'a- (2.35) 

The spatial components of the perturbed Einstein equations, imply, instead 

tP" + n{<P' + 2^') + {2n' + n'')<P + \v\<P - ^) = ATxGa\5,p, + 4pb), (2.36) 

V\<P -^Ij) = l2TxGa%p, + p,)VV, + {p, + p^) Wb], 5sPb = (2.37) 

In Eq. f l2.37p V^o"^ is the neutrino anisotropic stress, while V^o"b is the magnetic field anisotropic 
stress defined as: 

where, VL-q{x) is the magnetic energy density referred to the photon energy density and it is constant 
to a very good approximation if the magnetic flux is frozen into the plasma element. In closing this 
section it is appropriate to mention that, in the past, different approaches to the effects of large-scale 
magnetic fields have been proposed. In this respect the interested reader can usefully consult the the 
review articles reported in Refs. [13], [H] and [15] . 
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3 Autocorrelation of the Faraday rate 



The autocorrelation of the Faraday rotation rate can now be computed by taking into account the 
finite thickness of the last scattering surface. This analysis will be performed both analytically and 
numerically. The present results complement and extend the analysis of Ref.[37j which differs from 
other attempts [l6l HTJ HH SHI [50l EH [52] insofar as the contribution of the magnetic fields is included 
at every step of the Einstein-Boltzmann hierarchy, from the initial conditions to the relevant dynamical 
equations. 

Faraday effect has been studied in the context of a uniform field approximation [16], HTJ HH HHl [50] 
which breaks explicitly spatial isotropy. As discussed in the previous section, it is not consistent to 
keep the magnetic field uniform and, still, to assign the angular power spectra as if spatial isotropy 
was unbroken. On the contrary it has been correctly argued that the magnetic fields should not be 
uniform but rather stochastically distributed over different comoving scales as in Eqs. (I2.ip - fl2.2l) 
[33] [37] [STl [52] . The physical situation is such that magnetic fields with comoving inhomogeneity scale 
of the order of the Mpc are uniform over the scale of the Larmor radius but are not uniform on the 
scale of the Hubble radius. 

In [51] the interesting observation has been that one could indeed compute, in the sudden decou- 
pling, the correlation function of Faraday rotation. A similar observation has been later made in [52] 
within the same approximations. One limitation of is that, when computing the Faraday rotation 
angle, the visibility function has been taken to be infinitely thin and the decoupling sudden. The latter 
approximation is known to be inaccurate at small scales where the Faraday-induced B-mode autocor- 
relations increase. One of the themes of the present section will be to relax the sudden decoupling 
approximation in the calculation of the Faraday rate. 



3.1 Sudden and smooth decoupling approximations 

The sudden decoupling approximation offers a reasonable estimate of the angular power spectrum of 
Faraday rotation but only for sufficiently large angular scales (i.e. for sufficiently low harmonics since 
6 ~ vr/£). Conversely the smooth decoupling approximation is relevant for large i (in practice for 
i 3> 150) where the finite thickness of the last scattering surface (together with the other sources of 
diffusive damping) determines the fall-off of the Faraday autocorrelations. In the sudden approximation 
the visibility function is sharply peaked around the recombination time with zero width, i.e. /C(t) = 
6{t — Tree)- The rotation angle of Eq. fl2.30p can be expanded in ordinary spherical harmonics, i.e. 

Hn, To) = Y: fimYeUn), = E {fL fim), (3.1) 

£m m=—£ 

where C^^^^ denotes the autocorrelation of the Faraday rotation angle in full analogy with the au- 
tocorrelations of the polarization observables introduced in Eqs. (ll.lOp - Ql.lip . By assuming that 
/C(r) = 6{t — Tree), Eq. (12.301) implies that 

($(n,ro)<l>(m,ro)) = ^-^Cf^ Piin ■ m), (3.2) 



Cr = A\v)As{k^ToY--n'/' , ^ , / ^ (3-3) 
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Note that Eq. fl3.3p is derived under the approximation that 2{i + 1) < 5 — < (see [59] p. 675). 

When the finite thickness effects of the last scattering surface are taken into account the visibihty 
function can be approximated by a Gaussian profile centered at Tree, i.e. 

/C(r) = Ar((Trec)e / /C(r)t^r = l, (3.4) 



71 a. 



rec 



T0_ 



(3.5) 



eifiz) = ^ / e-'^dt. (3.6) 



a/tt Jo 

The overall normalization A/'(crrec) has been chosen in such a way that the integral of /C(r) is normalized 
to 1: the visibility function is nothing but the probability that a photon last scatters between r and 
T + dr. Equation fl3.5p simplifies when Tq ^ t^^c and Tq ^ (Xj-ec, since, in this limit, the error functions 



go to a constant and /^{crrcc) <^^cc y'^/'^- ^^e latter limit, the thickness of the last scattering 
surface, i.e. cXrec, is of the order of Tree- In the explicit calculations the normalization constant is always 
taken in its full form (see Eq. (13. 5p ): however, to avoid lengthy expressions the final results will be 
expressed, when feasible, in the limit tq ^ Tj-ec and tq ^ (Tree- 

Recalling that, in Eq. (I2.30p the magnetic field appears in real space, the $(n,ro) can also be 
written as: 

'^(^' ^o) = / ^'^^(^) ^(^) e-'"'^^'^-^^ dr, (3.7) 

where the present position of the observer coincides with the origin and where n = k ■ h denotes, as 
usual, the projection of the momentum of the photon over the Fourier mode. Inserting Eq. (13.51) into 
Eq. (13. 7p and performing the required integration over time the angular power spectrum becomes 

Cr = 4vr^'(^) + 1) / f ^B(fc)^^|^^'(A;, Tree), lC{k, r^ec) = e-'^^/'^S (3.8) 
where je{z) are the spherical Bessel functions [60], i.e. 

= ^J^Je+i/2iz). (3.9) 



To derive Eq. (13. 8p it has been used that, thanks to spatial isotropyo, {k ■ nf = A;^/3 and where we 
defined k^ = VQ/ctcc- The result of Eq. (13. Sp follows directly from Eq. (13. ip . More specifically, it can 
be shown that the fim are determined to be 

hm = ^'[2^)tt' -^(^)\/^(^+l) / d'k^-j^m . Y}^*Cmk, Tree), (3.10) 



where Y^rnW ~ ^fc^^m(^) / \J + 1) is the vector harmonic which arises from the analog of the 
Rayleigh expansion but in the case of a (solenoidal) vector field [HH [67] (like the magnetic field). 
Inserting Eq. (I3.10p inside Eq. (13. ip . Eq. (13. 8p is swiftly recovered. The angular power spectrum can 
then be written, in a compact notation, as 

^(FF) ^ ^(FF) ^ ^^^^^^ ^^^^ ^g^^^^ 

.-.4..,.,..(^)"-'..0.x.0-^^(-|)(^)-(^.)»-\ ,3.., 



^°Recall that, by assumption, the stochastic background of rehc magnetic fields does not break spatial isotropy and, 
therefore, n^n^ — (5*-'/3. 
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The integral appearing in Eq. (13. lip is defined as 

UnB,K,kn,h) = rx^«-^j2(xo)e-2(-«/^<^)' dxo, (3.13) 

JO 

where, on top of kt (aheady encountered in Eq. (13.81) ) kj) and k^ parametrize, respectively, the effects 
of the thermal and magnetic diffusivities [HE], i.e. Xd = k^TQ where 

kf = k^' + kf + k-'. (3.14) 

Note that, in Eq. (I3.12p . the normalization is fixed by choosing, as pivot frequency, z/max which 
maximizes the spectral energy density of the CMB (see Eqs. (I2.27l) - fl2.28l) and discussion therein). 
The three separate contributions arising in Eq. (I3.14p will now be estimated. 

As discussed in Eq. f l2.32p the maximum of the visibility function is determined by the condition 
e" + e'^ ^ which implies, in terms of the ionization fraction, that 

— 2HXe + fieaThxl ~ 0. (3.15) 

Equation fl3.15p defines Tree, indeed, at Tree, x'^^^ ~ — o"Th«(Trec)?^e(Trec)- The ionization fraction can be 
independently obtained from the approximate solution of the equation for recombination [5i 



.(..^i^(|if)"^(^)-'e-, ...3x:0^. ,3..) 

A slightly different parametrization can be found in Ref. [55]. From Eqs. fl3.15p - fl3.16p it follows that 
the visibility function peaks, within the present approximations, for z,.ec — 1050 and that Xg(ri.ec) — 
7Yrcc(-B/2;rcc)o"Th^c (Tree )fl (Tree)- Introducing therefore the parameter at — \/3B / z^ec we will have that 
£t = ktTQ can be written as 



2cKt \/Zxcc I , 2^eq 



Xa 

where 



1 + ^, (3.17) 



Since we will always be assuming the fiat ACDM paradigm, Eq. (13.181) depends effectively only upon 
f^Mo- The numerical result for the integral can be written as Xa = ^^mo^ and 71 = —0.0858. 

The shear viscosity term (to a given order in the tight- coupling expansion) allows for the estimate 
of ki). To zeroth order in the tight-coupling expansion the monopole and the dipole of the radiation 
field obey, respectively, the following pair of equations 



A/n H n H k^ — n + —, tAio 

^° i?b + 1 15 Rb + l ^° 3(i?b + 1) 

k^ r ^li 
^3(i?b + l) 



, // , ^-Rb , , k"^ 



B 

O'B 7- 



(3.19) 



where Axh = In the absence of Axh, Eqs. (13.190 and fl3.20p can be derived from Eqs. fl2.18p and 
f l2.2ip by recalling that, by definition, the monopole and the dipole of the intensity of the radiation 
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field are related, respectively, to the density contrast and to baryon-photon velocity, i.e. 5^ = 4Aio 
and ^-yb = SkAii. Introducing the photon-baryon sound speed 

' 3[l + i?b(r)]' ^ 4p^ VO.023 A^rcc + l/ ^ ' 

in Eq. fl3.19p the variable can be changed from r to = c^^dr. The resulting equation is solvable 
the WKB approximation [31]. Since we are dealing here with the polarization observables we are 
specifically interested in the evolution of the dipole. In the most simplistic case, i.e. when only the 
adiabatic mode is present 



An(fc,r) = _ v^^*W sin[fcr3(r)]e-^^/^B, r,(r) = T dr' c,^{t'), (3.22) 
5 JO 

where TZ^,{k) is a stochastic variable obeying Eq. (12.31) and where /cd can be estimated as: 

1 2 r , dr' 



.2 - , / Csb(r')^^^ . (3.23) 

/££) 5 Jo Ct(T j a^c^^cfTh 

The estimate of Eq. fl3.23p can be further improved by going to second order in the tight-coupling 
expansion where the inclusion of the polarization fluctuations allows for a slightly more accurate 
estimate 



dr' 



16 ^ B 



15 i?b + 1 



(3.24) 



kl Jo 6(i?b + l)e' 

The factor 16/15 arises since the polarization fluctuations are taken consistently into account in the 
derivation. To obtain explicit estimates of the diffusive damping it is very practical to use an explicit 
form of the scale factor. It can be shown by direct substitution that, across the matter radiation 
equality the solution of Eqs. fl2.5p - (l2.6l) is given by a(r) = aeq[(T/Ti)^ + 2(r/ri)]. To zeroth-order in 
the tight- coupling expansion, after a lengthy but straightforward calculation, we will have 

1 _ 9 gg ^ ^Q-sfhln,oY'^^fhlnuoV/'fio5iY/' 



knT, 



D ' roc 



0.023/ 



This estimate can be obtained from Eq. fl3.23p by setting i?b = in the baryon-photon sound speed. 
More realistically we have instead that 



Rdr) = i?b 



r \ ^ r 
- +2- 

nJ Ti 



:Rb = 0.216 f^V^^) . (3.26) 



In this case the estimate changes quantitatively but it can still be expressed as in Eq. f l3.25p but 
with a different pre-factor, i.e. 9.63 x 10^'^ becomes 2.17 x 10~^. More accurate expressions can be 
obtained from the second-order treatment by performing numerically the various integrations and by 
using, repeatedly, the interpolating form of the scale factor. 

The further source of damping contemplated by Eq. f l3.14p is the one connected with the finite value 
of the conductivity. In the electron-ion plasma characterizing our physical system the conductivity is 
determined by Coulomb scattering, i.e. 

3 / \ 1/2 1 

acoui = ^lnAc, Ac = - — -. (3.27) 
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To compute the conductivity we can assume that each colhsion produces a large deflection. Under 
this approximation we will have that 

T f T \ 

o= . =a,a, a, = — , (3.28) 

4 7rrcoui aemlnAcVme/ 

where Fcoui = "^e t'thO'coui is the Coulomb rate. The typical conductivity scale is then determined 
according to 

and it is typically smaller than the diffusive damping. The Silk damping is often introduced by 
summing up the two dissipative contributions, i.e. by defining £g ^ = i^'^ + which also implies that 

n 

. (3.30) 



3.2 Explicit forms of the Faraday autocorrelations at small angular scales 

The integral of Eq. fl3.13p can be evaluated in a closed form in terms of generalized hypergeometric 
functions 159): 



"B 3 I 
2 2 



x,(nB, fed) = 2-3^-^/2 {k^T.r^---^^ — ^ 



ai, a2, 03; ^1, ^2, 7. — 



(3.31) 



where, by definition of generalized hypergeometric, two of the three arguments are replaced by gener- 
alized vectors of dimensions q and p, i.e. Fgp{d; b; z) = Fqp{ai, ap] 61, 6^; z). In the case q = p = 3 
and for the case of Eq. (13.311) the components of the two generalized vectors are given by: 

a,=i+l, a2=i + ^, a3 = £+y-| (3.32) 

bi=i + l, ^2 = ^ + 1 63 = 2^ + 2. (3.33) 

For the present purposes Eq. f l3.3ip is not so revealing. More explicit expressions can be obtained by 
working out the large i limit of the spherical functions ji{xQ). In the limit i ^ 1 the spherical Bessel 
functions can be expressed as [6OI [61] 



fei^w) ^ r'm^T^l ^ '^(«^' ^) = ^v/^^^ - ^ arccos {w-') - (3.34) 
where w > 1 and Xq = iw. Using Eq. fl3.34p inside Eq. (13.131) it is easy to obtain: 

Ji \/vfi - 1 

Changing integration variable as w = a/?/^ + 1 we do get that 

/■oo — 2-^('w^+l') 

X,(nB, 4) = r^-^ / dy {y^ + l)("B-6)/2g ij^ >^ (g gg^ 

JO 
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where the oscillatory contribution is, in practice, averaged so that it is fully legitimate, in the large ^ 
limit, to demand cos^[/5(w, £)] 1/2. The integral over y arising in Eq. fl3.36p can be performed in 
explicit terms and it can be always expressed as 



where 



TC- 



5-nB 
2 



6-nB 
2 



1 nB_ 
2' ^ 



{2z 



,(5-nB)/2_ 



WB-5 
2 



r 



6— nt 



'11 



6 — ub 7 



:^,2.^ 



(3.37) 



(3.38) 



In the last equation, -Fii(c, d, z) denotes the confluent hypergeometric function [601 [61]. As a function 
of z, Q{z, ub) is slowly varying and it is always of order 1 for a variation of z between, say 0.001 and, 
e.g. 10 which is the physical range of harmonics. The result expressed by Eq. (13.381) seems to be 
naively singular when ub is a natural number since, in this case, the Euler Gamma function has a 
negative (integer) argument. This divergence is just apparent since by going back to the integral of 
Eq. (13.361) it can be explicitly checked that the result is convergent (both in the origin and at infinity) 
also when hb is a natural number. As an example consider the case ne = 1; in this case the result 
of the direct integration appearing in Eq. (13.361) can always be recast in the form of Eq. (I3.37P but 
where 

g{z, 1) = -zV'[2z2iro(^') + (1 - 2z')K,{z% (3.39) 

where Kq{z) and Ki{z) are analytically continued Bessel functions. In spite of the presence of the 
(increasing) exponential factor, Q{z,l) approaches zero for large z. For 2 — 0, Q{z,l) 2/3. It 
is amusing to notice that this result is consistent with Eq. (I3.38p . Indeed, in the limit 2; — >■ 0, 
Q{z, Hb) = v^/2[r((5 — nB)/2]/r(3 — nB/2) which equals 2/3 for ub = 1. 



4 Semi-analytic estimates of the polarization observables 

As discussed in the introduction, the calculation of the B-mode polarization is done, for our semi- 
analytic purposes, by iteration. In the first part of the present section, by applying previous methods 
described in [331 El], we ought to derive a semi-analytic expression for the E-mode autocorrelation 
which could be extrapolated for large multipoles with reasonable accuracy. In the second part of this 
section, the results of Section [3] will be combined with the estimates of the E-mode autocorrelations 
and we shall be able of obtaining a closed expression, valid at small angular scales, for the B-mode 
auto correlations . 



4.1 E-mode autocorrelations 

From the known expressions of the polarization observables (see, e.g. Eqs. (I1.7p - (ll.l0p ) it is easy 
to see directly that the presence of a B-mode polarization is directly related to a mismatch between 
A_|_ (n,r) and A_(n, r). This statement is clear in real space where 

AE(n,r) = -i{Al^)(n)[Ki')(n)A+(n,r)] + Kl-^)(n)[iri-')(n)A4n,r)]}, (4.1) 
AB(n,r) = ^{iri^)(n)[i^i')(n)A+(n,r)] -;^l-^)(n)[iri-')(n)A4ri,r)]}. (4.2) 
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The differential operators appearing in Eqs. (14. ip and (14.21) are generalized ladder operators [H [7] 
whose action either raises or lowers the spin weight of a given fluctuation. They are defined, in general 
terms, as acting on a fluctuation of spin weight s: 

1 



Klin) 
K' (n) 



1 

(sin-i?)*^ 



smv 



[sin 

(sin^9)^ 



(4.3) 
(4.4) 



For instance, K^^A^ will transform, overall, as a function of spin-weig ht 1 while K^y[K^yA+] IS, as 
anticipated, as scalar. Using Eqs. (14.31) and (14.41) inside Eqs. (14. ip and (14. 2 p the exphcit expressions 
of the E-mode and of the B-mode are, in real space: 



AE(ri,r) 



dUA^ + A_ 



+ 2i 



d^d.iA, 



) -2(A+ + A_) 
,9^(A+-A_)]}, 

AB(n,r) = ^{(l-^2)^J(A+-A_)-4/x9^(A+-A_)-2(A+-A_) 



1 — fi- 



I 
2 

dl{A, 



A. 



+ 2i 



d^d^{A+ + A_ 



r5^(A+ + A_ 



(4.5) 



(4.6) 



In Eqs. (14.50 and (14.60 9^ denotes a derivation with respect to /i = cos while (9^ denotes a derivation 
with respect to (p. Equations (14.51) and (14. 6 p hold in the spherical basis where n = (t?, v^). From Eqs. 
(14.50 and (14.60 it can be concluded that the B-mode polarization is absent (i.e. Asin, t) = 0) provided 
A+(n, r) = A_(n, r) and provided i9<^(A+ + A_) = 0. The latter condition implies that the brightness 
perturbation is independent upon cp. Indeed, to lowest order in our iterative scheme we will have that 
A+(n, r) = A_(n, r) = Ap(n, r). Moreover, in Fourier space Ap(A;, /i, r) does not depend upon ip and 
it obeys 

= (Ai2 + Apo + Aps), (4.7) 



A'p + {tkfi + e')Ap = —(1 - fi^)Sp, 

where Sp{k,T) is computed consistently with the other evolution equations which include the large- 
scale magnetic fields both at the level of the initial conditions as well as at the level of the dynamical 
equations. To lowest order, the heat transfer equation does not depend upon p the following pair of 
equations will hold, in real space: 



AE(n,r) = -d^[{l - iJ, )Ap(n,r)], a2, em = a-2, £m- 
As a consequence of the second relation in Eq. (14.81) . it also follows that 



,(E) 
''em 



dn Ap{x, n, ro)[2F;„(n) +_2 Y;^{n)], 



smv 



2? 1 

sm v ^. 



Y, 



tm 



n] 



(4.8) 

(4.9) 
(4.10) 



where ±Yem{n) are the celebrated spin-2 spherical harmonics. Using line-of-sight integration, the 
autocorrelation of the E-mode polarization (defined in Eq. (II. lip ) can be written as 

dk , 



^(EE) 

^e 
Af\k) 



4- / ^-^Vn{k)\Af\k)\^ 
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(4.11) 
(4.12) 



The E-mode autocorrelations of Eq. fl4.12p can be estimated with semi-analytic methods. For this 
purpose the first step is to derive Sp{k,T) in the tight-coupling approximation. This first step is 
known to be insufficient for a reliable estimate of the E-mode autocorrelation. Therefore, following the 
approach originally developed in [53] , the tight-coupling expansion can be systematically improved by 
deriving an effective evolution equation for Sp, i.e. 



S^ + ^e'Sp = klAn{k,T) 
10 5 



(4.13) 



where Au{k, r) is the dipole of the intensity (related to the velocity of the photons) and evaluated to 
zeroth-order in the tight-coupling expansion. Using the solution of Eq. fl4.13p we then have 



Jo Jo Ji y 



dy 

y 



Recalling Eq. f l3.22p we can express the angular power spectrum of the E-mode polarization as 



(EE) 



C 



+ !)(£ + 2)(£ - 1) {j-y {k,a,,,fX,{n, 



(4.14) 

(4.15) 
(4.16) 



where C = 3.5 x 10^'^ ^t? and where 



7(^rec) 



1 



To JO 



V3to Jo 



Ticc/ti 



dx 



1 + /?(a;2 + 2x) 



(4.17) 



The integral appearing in Eq. ( 14.16^ can also be written as: 



s, 



(9ns -5 



2 
1 



(^2^1)K-6)/2g-2(^V^^)(,/+l) 



{Us, £d) = ^r=-^ / {y' + 1)("=~^)/^ COS [2£7(r.ec) Vy^]e-^(^^/^o)(.^+i) 



It turns out that a good approximation to the full result can be obtained by setting 

ris-l 



c, 



(EE) 



^(EE) f « 



(4.18) 
(4.19) 

(4.20) 
(4.21) 



where £p is a given putative frequency. The estimate of the angular power spectrum of the E-mode 
can be improved by solving more accurately Eq. f l3.20p which can be written as 



Aii(A;, r) = - cl!^[aTi{k, r) - a^{k)] sin [kr^{T)\. 



(4.22) 



where 



an{k,T) = ^T^*{k) In 



25^'"^''. 



aBik,T) = 3^/^ 



ik)-Mk) 



R 



+ ^^B{k) \e 



(4.23) 
(4.24) 



These expressions allow for more accurate results. However because of the logarithmic factors they 
are more difficult to handle inside the integrals. 
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4.2 B-mode autocorrelations 



Having found the E-mode polarization to zeroth order in the Faraday rotation rate, the calculation 
can now be extended to first-order, i.e. when A+ 7^ A_. As previously noticed [37], for magnetic 
fields of comoving intensity smaller than 10 nG the Faraday autocorrelation can be safely evaluated in 
the limit of small rotation ang ifl This will also be the strategy followed here so that, in real space, 
the first-order iteration leads to an expected mismatch between the two aforementioned orthogonal 
combinations: 

A+(ri,r) = -2<l>(ri,r) Ap(f,n,r), A_(fi, r) = 2 <l>(n, r) Ap(fi, r), (4.25) 

where Ap(n, r) follows from the solution of the heat transfer equation to zeroth order in the Faraday 
rate. From Eqs. (14.51) and (14.61) the B-mode polarization turns out to be given by: 

Ab(x, n,r) = 2dl[il- fi^)^x,h,r)Apih,r)] (4.26) 

«S = ^/rfnF4(n)AB(n,r), (4.27) 
whose explicit expression can also be written as 

aS = Jdn $(n) Ap(ri) [,Y;^{n) Y;jn)]. (4.28) 

For sufficiently small angular separations, the 2-sphere degenerates into a plane where the fiat basis 
is characterized by a pair of (Cartesian) coordinates and "^y. In the fiat basis (as opposed to 
the spherical basis) the third direction is fixed. The limit in which the 2-sphere degenerates into a 
plane is one of the simplest examples of group contraction [58]. Owing to this observation the matrix 
elements of the representations of 5*0(3) can be related to the matrix elements of the representations of 
M(2) (i.e., in the language of |5^, the 5*0(2) group supplemented by the two-dimensional Euclidean 
translations). In the limit of large multipoles, therefore, the Legendre polynomials degenerate into 
Bessel functions, i.e., more precisely [58] : 

hm p/cos = Mr) ~ MM). (4.29) 

At the right hand side Jo{r) is a Bessel function of zeroth order while at the left-hand side it appears 
a Legendre polynomial whose argument can also be written as P^(cos't?). The limit (14.291) corresponds 
to keep fixed M by sending i ^ 00 and —* 0. 

Every calculation conducted in the spherical basis can be performed in the fiat basis and the 
results of the two calculations must coincide in the limit i 00 and 6' — > 0. Consider, for instance, 
the temperature fiuctuations computed in the spherical basis 

Ai(#) = 7^ / d^£Ai{£) e''-\ (Ai(f)Ai(f')) = C^i^^^6^^\I+ £'). (4.30) 
A simple calculation in the fiat basis shows that 

OTT(t9i2) = (Ai(/i) Ai(/i)) = ^ hM CrMOui), (4.31) 

^^The rationale of this statement is very simple. If we consider the smallest imaginable frequency we can derive which 
is the bound on the comoving magnetic field over the pivot scale corresponding to a Faraday |$(?i)| < 1- The result 
of this analysis [37] imply, roughly, that for comoving frequencies of 30 GHz, comoving fields of 100 to 1000 nG still 
lead to a rotation rate smaller than 1. For larger frequencies the bound becomes looser. For a range of variation of riB 
between, say, 1 and 2.5 the constraint changes by less than one order of magnitude. 
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where the Bessel function arises from the integral 



47r2 J ' ' J ~ Jo '^"^ ' 27r 

The second equahty in Eq. ( 14.32^ follows from a well known integral representation of Bessel functions 
[HDl EI]; the result of Eq. 04.321) can also be obtained from the temperature autocorrelation computed 
in the spherical basis by consistently taking the limit i oo and ^9 — > 0: 

c^^{A2) =j:^^cr^pd^osA2) - 7^ / idicrMe^^i), (4.33) 

where the second expression follows from Eq. f l4.29p . Equations f l4.30p -( l4.33p demonstrate the equiva- 
lence of the two computational schemes. The polarization observables and the spin weighted spherical 
harmonics can be expressed in the flat basis [HI [7] (see also [621 [63] for a connection with the well 
known problem of computing ellipticity correlations from weak lensing in the fiat basis). 

For a semi-analytic estimate of the B-mode autocorrelation, Eq. (14.281) can be inserted inside Eq. 
(11. lip . The derivation will be performed in the spherical basis, mentioning, when relevant, how the 
results also emerge by working in the fiat basis. Following the calculation originally due to Goldberg 
and Spergel [61] the B-mode autocorrelation can be recast in the form 



^ jdj, /rf72-^(£,Ji,j2)C7jf)(^D,A)C7jf^^(^D,A), (4.34) 

where sums have been replaced with integrals since, in the small-scale limit, ji and j2 become, in 
practice, continuous variables. In Eq. (14.340 the dependence of the angular power spectra upon 4 and 
£d has been made explicit. The exact form of JF(£, ji, J2) can be written as: 

^1^,J1,J2J 4^ (,2£+l^(£-l)(£+l)(£ + 2)j2(j2-l)(j2 + l)(j2 + 2) ' ^^-"^^ 

where 



/(£,J1,J2) = 2jij2(jl + l)(j2 + l)+2ji£(ji + l)(^+l)- £2(^+1)2 

- j\2(ji + If - j|(j2 + 1)' + 2j2(j2 + 1) + 2£(£ + 1) - 2ji(ji + 1), (4.36) 

Kono? = -^j_^d^^Peifi)P,,{fx)P,,{fx). (4.37) 

In quantum mechanics, the Clebsch- Gordon coefficients are probability amplitudes arising when we 
have two different angular momenta (be them Ji and J2) acting on two different Hilbert spaces. The 
complete and orthonormal set of eigenvectors of and of Jiz are |ji mi). Similarly, the complete and 
orthonormal set of eigenvectors of J| and of J22 are 172^2). By doing the tensor product between 
the vectors of the two subspaces with (2ji + 1) and (2^2 + 1) dimensions we obtain the decoupled 
representation whose basis vectors are characterized by four different eigenvalues, i.e. ji, j2, mi 
and 7712- The same problem can be addressed by exploiting the coupled representation where the 
eigenvectors are characterized by a different system of eigenvalues i.e. i, m together with ji and j2- 
Since both systems of eigenvectors are complete and orthonormal they are connected by a unitary 
transformation 

\jlj2im) = (^ii mi j2 m2 lil J2 mi m2 ) , 



mi, m2 

2 

-'jl mi J2 m2 I 



IC'i'^^i ,2 m2r = I0U2 mimsbi J2 ^m)\\ |ji - J2I < ^ < Ji + J2, (4.38) 
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reminding that the Clebsch- Gordon coefficients vanish unless the triangle inequality is obeyed. From 

■^l m 



Eq. fl4.38p it can be easily argued that ^2 Wigner symbol up to a numerical factor 



(i.e. V^e + 1) and up to a phase (i.e. (-l)^i-^2+™3)[65]. Furthermore, from Eq. (gSHD, Cj^ojaO is 
zero if the sum of the eigenvalues gives an odd integer and and it is different from zero otherwise. All 
the mentioned properties can also deduced in the particular case of Eq. fl4.37p and bearing in mind 
the recurrence relations of the Legendre polynomials. The limit where ji ^ 1, j2 ^ 1 and i ^ 1 
is physically opposite to the "quantum" limit where Eq. (14.381) originally arises. We could say that, 
using the quantum mechanical analogy, the small-scale limit is "classical" since there the Clebsch- 
Gordon coefficients can be thought to arise from the composition of classical angular momenta. This 
approach to the asymptotics of the Clebsch-Gordon coefficients was originally studied by Ponzano and 
Regge [68] by exploiting the connection of the Clebsch-Gordon coefficients with the Wigner 3j and 6j 
symbols. In the case which is mostly relevant for the present investigation we have that [68] (see also 
[691 [70]) 

IC^O |2_ ^^+^ ^ (Ann) 

It is worth mentioning that, in the general case, the semiclassical limit of the Clebsch-Gordon coeffi- 
cients [69] (see also [70]) has been studied with a technique which is closely related to the perspective 
adopted in the present investigation. Indeed in [69] the results of [68] have been con&med 1^ by using 
the recursion relations obeyed by the 3j and 6j symbols. The latter relations become, in the semi- 
classical limit, differential equations. In the present approach we somehow use the same observation 
when replacing discrete sums by integrals (see Eq. fl4.34p ). In the flat basis, the result of Eq. (14.390 
follows from the evaluation of the integral 

POO 

n Mhy)Mhy)Miy)ydy. (4.40) 



which coincides exactly with the right-hand side of Eq. (I4.39P (see, e.g. [581 [59]) when \ii — £2! < ^ < 
ii + £2 and it is zero otherwise. Inserting Eqs. (I4.36p . (I4.37P and (I4.39P into Eq. (14.350 we do get, 
after simple but lengthy algebra, that 

(4.42) 



(x + y) [l 


^2 (x + yf-2x^ y2_2£4 


2 


27r2 £4 (a; _ yf 

(x + y) 


^{x - £) (a; + £) (£ - ? 
' {x + y f -2S:^f -2' 


2 


27r2 (x-y)^^ 


/{x-1) (x + l) (l-y) (l + y)' 



In Eq. (I4.4ip . for immediate convenience, the following change of variables has been adopted 

x = ji+ J2, y = 31-32, 0<y<i, x> i, (4.43) 

where the range of variation of x and y follows from the triangle inequality. In Eq. (I4.42p a further 
change of variable has been performed, i.e. x —>■ x = x/i and y ^ y = y /£. In terms of x and y the 
integration region becomes < y < 1 and x > 1. 



^^It must be said that the discussion of [68] was, in its nature, rather formal. The authors defined asymptotic 
expressions vahd in quantum, classical and semi-classical domains without discussing the smooth connections of the 
different limits. This aspect has been later deepened in [Ml ED] . This is why it will be important to cross-check the 
obtained results with direct numerical evaluations of the relevant matrix elements. This will be one of the themes of 
the following section. 
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The largest signal obtainable from Faraday rotation arises for frequencies smaller than the maxi- 
mum of the CMB (i.e., roughly, v < i/max — 222 GHz) and for angular scales much smaller than one 
degree. Over these scales C^^^^ increases with i. Recent experiments Quad [15], [16], [17] give upper 
limits on the B-mode polarization also for i ^ 200 and, with large error bars, up to £ ~ 1987. The 
autocorrelations of the polarization and of the Faraday rotation can be expressed, respectively, as 

Cr = A/'f nB)e-^^^/^^ = ATe r^-'g{i/in, n.)e-'''/'i, (4.44) 

in the small scale limit. The results of Eqs. (14.431) and (I4.44p then allow for the following estimate of 
the B-mode autocorrelation 

«*+"ci-'.^BriV™"'-""%-»'''^->' (4.45) 
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where ai = 1.13 and where 



= £^' + if^2£^' + £;' + if, (4.46) 
ATb = ^ATpAr^C'"" ) ' ( ) ^ ( , (4.47) 



-6 /VAb\V ^ \-V^0^"^"' 



AT. ^ 5.075 X 10- ( ^ j ^^J , (4.48) 

AAe = 2.96 X 10-^ (^y' (/.K)l (4.49) 



,2.41 X 10-V Vfcp/ V2.725 
The result of Eq. fl4.45p follows from the evaluation of the integral: 

'W{u,v,nB,ns) = / dx dy{x + yY^'"^ {x - yY'''^ x 
Ji Jo 

X [{x + yf-2x^f-2]^ ^~u^(i+y)y2 ^-vH,-yf/2 (4 
^(5-l)(£ + l)(l-y)(l + y) 

where u = i/id and v = i/i-o- As indicated >V(u, v, ub, Us) depends upon four parameters, i.e. tt-b, ''^s, 
u and V. Since we will be mainly interested in the case where rig is closed to its best-fit value we shall 
assume that rig varies between 0.9 and 1.1. The allowed variation on n-Q must necessarily be broader 
and it will be assumed, quite generally, that < ne < 5 which include the case of red (i.e. < I), 
scale-invariant (i.e. ne — 1) and blue (i.e. ne > 1) magnetic spectra. Finally one can also assume that 
u and V are of the same oder, i.e. they can be different but their mutual deviations are well within the 
order of magnitude and, typically, even less than a factor of 2. In this restricted parameter space the 
full integral of Eq. (14.501) can be reasonably approximated by the following expression: 

W{U, V, riB, n,) ^ 2"s+"=+=^/2(^2 ^ ^2y/2-n,/4-ns/i ^.a^iu^W)^ IQ~3 < u,v < 10. (4.51) 

Note that in Eq. (14.511) the range of variation of m ~ f is rather broad and much larger than what is 
needed. Equation (I4.5ip is certainly rather handy and shall be used, in a moment in direct estimate. 
However, an even better analytical approximation can be obtained by integrating directly Eq. (14.501) 
for a given set of parameters. Before switching to the phenomenological implications of the present 
results, it is relevant to remind that the present estimates are all conducted in the framework of the 
(minimal) ACDM scenario, in this sense the restrictions on the variation of the scalar spectral index 
are fully justified. 
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5 Parameter space and parameter extraction 



The expression of Eq. fl4.45p attains its maximum for 



1 Us + riB + 2 - 

2 V «i 



(5.1) 



In Eqs. flOTD and fOSD . following the WMAP5-yr data alone, ko = 7.08 x 10"^ Mpc"^ The value 
of the pivot scale kp is set, as usual to 0.002 Mpc~^. The value of the magnetic pivot scale /cl will be 
taken to be lMpc~^. The estimates of the previous sections will give, for J^, a value laying between 
i ~ 900 and i ~ 1100. In the case of the best fit parameters of the WMAP5 data alone we will have 
that Zd — 911. In Fig. [2] the full numerical calculation is compared with the analytic results obtained 



B, = 5 nG, n„ = 2, n =0.963, v = 1 00 GHz 




B, = 2 nG, rig = 1.5, =0.963, v = 100 GHz 




Figure 2: Comparison between the analytic and numerical estimates of the magnetized B-mode auto- 
correlation. 

in the previous section. In the plot at the left the magnetic field has been set to -Bl = 5 nG (with 
n-Q = 2) while in the plot at the right i?L = 2 nG (with ne = 1.5); i?L is the value of the magnetic 
energy density regularized over a spatial domain L ~ /^l^ where /cl = Mpc~^ is the magnetic pivot 
scale introduced in Eq. (12. 2p . The regularization scheme depends upon the spectral slope: for ne > 1 
(blue spectra) it is often performed by means of a Gaussian window function. In the case of red spectra 
(i.e. ub < 1) the regularization is done by means of a smooth representation of a step function. The 
relation between and Ab is given by 

Ab = ^ ^ Bl nB>l, (5.2) 



r( ^ 



1 — Ub f h 



1— riB 

Ab = ^^[^) Bl nB<l. (5.3) 

The passage from Ab (which has dimensions of nG^) to i?L can be done at any stage of a given 
calculation and it is not really crucial. For instance, in the case of curvature perturbations the fit to 
the experimental data is customarily performed in terms of An which is the analog of Ab and not of 
i?L- In both plots of Fig. [2] the full line is obtained from Eq. (14. 44 p . The dashed and dot-dashed lines 
are instead the result of a direct numerical integration using the code developed in [351 [37]. In 
the plot at the left the theoretical result has been divided by a factor 2 while in the plot at the right 
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the result of Eq. fl4.44p has been multiphed by a factor 2. This is due to the fact that the integral of 
Eq. fl4.49p has been approximated with Eq. (14.501) . In other words the approximation of Eq. (14.501) is 
rather good for the slope at the expenses of a factor of 2 in the amplitude. The latter statement can 
be better understood by noticing that the position of the maximum is always well determined by Eq. 
(15.11) and this reflects the correctness of the £ dependence in Eq. (14.441) . Equations (I4.49p and (I4.50p 
are quite useful if we wish to set direct bounds on the magnetic field intensity starting from published 
upper limits on the B-mode polarization. 

We are now going to compare the obtained results with the direct bounds on the B-mode autocor- 
relations. There are two complementary ways of making use of the data illustrated in Tab. [H The 
first strategy is direct, i.e. from the upper limits on the B-mode autocorrelations, a constraint on the 
properties of the magnetic field intensity can be inferred. The second strategy is to take the data of 
each experiment at face value and to perform a global analysis. Four of the experiments of Tab. [1] 
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Figure 3: Bounds on the parameter space of the magnetized background from Maxipol, Boomerang, 
Cbi and Capmap. The full line in the right plot denotes the foreseen signal for a stochastic background 
of relic gravitons with rx = 0.2. The shaded area denote allowed regions of the parameter space. 

reported direct upper limits on the B-mode autocorrelation. Maxipo l^p^ [H] and BoomeranJ^ [T2] 
derived, respectively, the following upper limits on the B-mode autocorrelation: 



^{^ + 1) 



Cf^) < 112.36 {^iK)\ 



+ 1] 



Cf^^ < 8.6 (/iK)' 



(5.4) 



where both figures are at 95% confidence level. The Maxipol limit becomes more tight (i.e. from 
112.36 (yU K)^ down to 90.25 (yU K)^) if the calibration uncertainty is excluded [131 [H]. Both Boomerang 
and Maxipol operate at a relatively high frequency (i.e. 140 GHz for Maxipol and 145 GHz for 
Boomerang). Cbllfl [18j and Capmap0|W| gave more stringent upper limits, i.e., respectively 

£{£ + !] 



+ 1 



Cf"^ < 3.76 (yuK)^ 



Cf"^ <4.8 {fxKf 



(5.5) 
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The upper limits of Eq. (15. 5p are more stringent, in absolute value, than the ones of Eq. (15. 4p . 
Furthermore, as discussed in the previous section, the magnetized B-mode autocorrelation scales with 



^•^Maxipol is an acronym for "Millimitcr Anisotropy experiment Imaging array" 

^''Boomerang is ana acronym for "Balloon Observations of Millimetric Extragalactic Radiation and Geophysics" 

^^Cbi is an acronym for "Cosmic Background Imager" 

^^Capmap is the acronym for "Cosmic Anisotropy Polarization Mapper". 
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the frequency. Therefore the signal is larger for lower frequencies. The Cbi and the Capmap experiment 
are both working at frequencies lower than Maxipol and Boomerang. More precisely Cbi operates 
between 26 and 36 GHz. In the case of Capmap twelve receivers operate, approximately, in the W 
band (i.e. for 84 GHz < u < 100 GHz) while the remaining four receivers operate in the so called 
Q band (i.e. for 35 GHz < u < 46 GHz). In Fig. [3] the bounds stemming from Eqs. (15.41) and 
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Upper iimit from CBi (95 % CL) 
Upper iimit from CAPMAP (95% CL) 
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Figure 4: The same bounds illustrated in Fig. [3] are here illustrated in terms of i?L and for blue spectral 
index. As in Fig. [3] the shaded areas denote allowed regions of the parameter space. As in Fig. [3] the 
region below the full (thick) line in the plot at the right is not excluded but it would correspond to 
rather optimistic sensitivities of the forthcoming experiments. 

(15.51) are compared with the signal of the magnetized B-mode. Concerning the exclusion plots of 
Fig. [3] few remarks are in order. Since the experiments gave specific upper limits, they must be 
applied, by definition, where the signal is larger and this occurs, in our case, for low frequencies and 
for i ~ £max- The common putative frequency has been chosen, for the purposes of Fig. [3l to be 30 
GHz. Both plots of Fig. [3] contain the same amount of informations. In the plot at the right the 
bottom line denotes the equivalent signal on the B-mode autocorrelation in the case of a stochastic 
background of gravitons with tensor to scalar ratio ~ 0.2. The latter figure is typically one of the 
goals of forthcoming experimental programs aimed at observing the B-modes arising from the tensor 
modes of the geometry. In this sense the shaded area allows to estimate the region of the parameter 
space which will be eventually accessible when the given experiment will improve their sensitivities 
by attaining progressively lower values in r^. Since some authors insist in setting limits only for the 
case of blue spectral indices, in Fig. H] the upper limits of the various experiments are expressed in 
the plane (ns, -Bl) but only in the case of > 1. The exactly scale-invariant case (i.e. ub = 1) 
would lead to a logarithmic divergence in the two-point function in real space. This is the same kind 
of logarithmic divergence encountered in the case of curvature perturbations in the case of the (exact) 
Harrison-Zeldovich limit (i.e. ng = 1). In the latter case we do have that 



(7^(x,ro)7^(y,ro)) 



dk 
T 



sin kr 

kr 



7^■ 



dk sin kr 
k kr 



(5.6) 



where the second equality follows from the first by recalling that Vn = Anik / kp)"""'^ and by setting 
ris = 1. As in the conventional ACDM paradigm the point rig = 1 is not artificially excluded from 
the analyses, for the very same reason the point ne = 1 is by no means less physical than the region 
nB > 1 and must be considered when scanning the parameter space of the magnetized ACDM model. 
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n =0.963, 1 = 2-6 n =0.963, 1= 2-6 




Figure 5: The WMAP upper limit applied to the case of the magnetized B-mode. 
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Figure 6: Confidence regions using WMAPSyr data pTl [22l [23l[2il [2^ 



Also the WMAP collaboration published an upper limit on the B-mode autocorrelations and we are 
now going to compare such a bound with the magnetized B-mode. The WMAP collaboration states 
an upper limit at 95%CL which can be phrased as follows: 

<Q_^5(^j^)2^ 2<£<6. (5.7) 



27r 

Since to derive the upper limit only the Q and V channels have been used the bounds should be imposed 
both at 41 and at 61 GHz knowing that, however, the most constraining frequency will always be the 
smallest. In Fig. [5] the WMAP upper limit is applied to the magnetized B-mode. The resulting 
constraint is extremely weak and it is consistent with published bounds [37]. So far the upper limits 
provided by the various experiments have been taken at face value and the consequent bounds on 
the parameter space of the large-scale magnetic field have been derived. It is also possible to use the 
data on the B-mode polarization in a dedicated analysis. The idea pursued here is rather simple 
and it is based on the observation that, nowadays, the differences between bayesian and Gaussian (i.e. 
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Figure 7: Confidence regions using Quad data. 



frequentist) approaches are becoming less relevant. For instance, in the joint two-dimensional 
marginalized contours at 68% and 95% CL have been reported for various combinations of parameters 
in the case of the WMAP data alone (see, e. g. Fig. 10 in Ref. [29]). The various marginalized 
contours are ellipses which have, approximately, the Gaussian dependence on the confidence level. 

In this situation it is legitimate to consider all the standard ACDM parameters as Gaussian variables 
whose mean is given by the best-fit value, for instance, of the WMAP 5-yr data alone: 



(fibo, fico, ^^A, ho, n„ r) = (0.0441, 0.214, 0.742, 0.719, 0.963, 0.087). 



(5.^ 



Following this frequentist approach the parameters of the magnetized background can be estimated 
by minimizing 



X (nB,v4B) 




(BB). 



(5.9) 



where the sum goes over all data points; C^^^\i) and a{i) are the observational data and their 
binned errors. Finally C^^^\i,nB, A^) are the values of the magnetized B-modes computed using the 
approach outlined in the previous sections. The {tibjAb) - parameter space is explored using a grid 
of 450^ points. In the following the contour plots for 68% and 95% confidence level are presented for 
different data sets. Since the instruments operate at different frequency channels, however, the data 
are usually those for combined frequency channels, it seems useful to assign effective frequencies to 
each of the experiments. Following [71] an effective frequency i/es can be defined by 



eflfj 



f{u)Siu)du, 



(5.10) 



where 5'(z/) is the frequency spectrum of the signal, and /(z/) is the frequency response profile of the 
detector. The polarization signal is determined by AB{n) which is proportional to i>~'^, so that the 
spectrum of the signal is given by 

S{u) = So(^) \ (5.11) 
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Figure 8: Confidence regions using the combined data sets of Capmap, Dasi, Cbi and BOS. 

Tliere are different options for tfie detector response function. For a fiat bandpass filter witli sliarp 
cut-offs at some frequencies z/^ and vi, tfie effective frequency is given by z/gfj = ^JVaV\,. Applying tfiis to 
tfie Dasi and Cbi experiments gives in botfi cases z/eg = 31 GHz. Assuming that the detector response 
function is modelled by a 5 function at the relevant frequencies it is found that the effective frequency 
is given by v^^ = ^^''^^^ , which in the case of QUaD results in z/gf/- = 118 GHz and the Q+V data 

of WMAP in Veff = 48 GHz. For Capmap combining the two detector response functions an effective 
frequency of z/ = 52 GHz is found. 

We are now ready to discuss the main results of this analysis. In Fig. E] results are presented for 

1 

WMAP 5. The minimum value of corresponds to = 12.89 nC and ub = 0.35 which corresponds 
to a magnetic field strength Bl = 332 nC. Furthermore, Xmm — 5-8 and there are 8 data points. In 

Fig. [3 the results are presented for the Quad data which consist of 23 data points. It is found that 

1 

the bestfit model has Xmin = 11-65 and corresponds to the magnetic field parameters A]^ = 88 nG, 
Ub = 1.49. This implies a magnetic field strength = 108 nG. By comparing Fig. [7] and [6] it seems 
that the Quad data lead to systematically larger values of the magnetic fields. At the same time, by 
admission of the Quad team, the BB angular power spectra might not be free of systematics. Indeed, 
the argument used by the Quad team seems to be, in short, the following (see [I6j). If Quad would 
have seen BB power at a level greater than in the ACDM, extensive investigations would have been 
required to check if the result was indeed free of systematics. Since the collaboration did not see 
any excess in the BB power, "further investigation of mixings systematics is arguably unnecessary" 
(verbatim from [16]). The caveat is that, in the present case, the putative signal is larger than in the 
ACDM case and, probably, further checks on the systematics may resolve the problem. 

The results for the combined data sets of Capmap, Dasi, Cbi and Boomerang B03 are presented 
in Fig. [HI In total there are 24 data points. The bestfit model is characterized by Xmm = 12.93 and 
the magnetic field parameters, A]^ = 0.8 nG and ub = 0.76. This results in a magnetic field strength 
Bl = 6.23 nG. Concerning the results expressed by Figs. El H and [H] few comments are in order. In all 
cases the relevant statistical indicator, for the present analysis, is the Ax^, i.e. ~ Xmm- We do not 



29 



judge here the compatibihty of a model with the data by performing the Pearson test. 



6 Concluding remarks 

In summary, the B-mode due to Faraday rotation has been calculated both semi-analytically and 
numerically. The parameter space of the magnetized background has been investigated by enforcing 
the upper limits on the B-mode polarization and also by performing a dedicated analysis. The two 
analyses led to compatible results. It has been shown that forthcoming experiments aiming at a direct 
detection of a B-mode polarization arising from a putative tensor power spectrum can also be able to 
cut through the parameter space of a magnetized background. 

An interesting byproduct of the present analysis has been an analytic scrutiny of the Faraday 
induced B-mode at small angular scales. In the small-scale limit the B-mode autocorrelations can 
be computed with reasonable accuracy by contracting (in the Inonii-Wigner sense) the relevant ma- 
trix elements of the rotation group. The latter results have been cross-checked by direct numerical 
calculations. 
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